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It is well-known that the universe was in a plasma state both before decoupling and 
^ ' after reionization. However, the conventional wisdom has been that the plasma effects 

fv^ ' are largely Debye-shielded and can thus be safely ignored when considering large scale 

QQ evolutions. Recently we showed that large scale structure formation in the universe may 

fvj , actually be suppressed by the plasma collective effect. Indeed, observational data indi- 

f^ ■ cate that the conventional theoretical formula tends to overestimate the matter power 

spectrum at scales k > l/iMpc~^. In this paper, we further develop our theory through 
a more thorough and general derivation of the Maxwell-Einstein-Boltzmann equation. In 
addition to baryon density perturbation post reionization, we apply this general formu- 
lation to investigate the possible plasma effect on CMB anisotropy. As expected, while 
the plasma effect does render an observable effect to the former, its impact on the latter 
is totally negligible. 
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1. Introduction 



Recently we showed that large scale structure formation in the universe would 
be suppressed by the plasma collective efFectM^. Indeed, observational data seem 
to indicate that the conventional theoretical formula tends to overestimate^ the 
matter power spectrum at scales k > l/iMpc~ . It is implied that the underlying 
theory may have to be modified. Moreover, our worw^ is performed in the regime 
of Newtonian limit and we expect to generalize the plasma collective effect in the 
general relativistic regime. 
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On the other hand, the cosmological perturbation theory(CPT) is not complete 
in the sense that the electromagnetic interaction is absent. Einstein-Boltzmann 
equations accounti ng fo r evolution of cosmological perturbations have been ex- 
plored for decade d ' I I Concerning anisotropics, the Sachs- Wolfe effeclP and Silk 
damping^ were predicted just after the d iscovery of CMB"' and the approach of 
tightly coupled limil!^ was perfecte d ^ ^ in 1990's to account for the CMB spec- 
trum. Neither these studies include EM fields even though tightly coupled limit 
approach combines photons and baryons as a cosmic plasma, nor have investiga- 
tions been performed on the likeliness of plasma effects in CMB study. In parallel, 
the effor ts on ma. tter inhomogeneities including study on the growth and transfer 
function j I I I comprise the evolution under only gravity. Later modificatiorM^ 
taking baryons into account neither addressed EM interactions, though it deri ved 
more accurate fitting formula to replace the pure cold dark matter BBKS fitting^. 
Though the CPT as well as its predictions and description of CMB and MPS do 
not take the EM interaction into account, we expect plasma collective effects for 
the universe was in plasma state before decoupling and after reionization. 

In the paper, we generalize the Einstein-Boltzmann(EB) equations in the CPT to 
our Maxwell-Einstein-Boltzmann(MEB) equations, which describe the perturbation 
evolution under both the gravity and EM interaction. With the inclusion of plasma 
effects we derive not only CMB ion-acoustic oscillations but also baryon ion-acoustic 
oscillations. The paper is organized as follow. In Sec. II, we justify that the plasma 
conditions apply in the universe during pre-decoupling and post-reionization era. 
Based on this justification, we generalize the EB equation to the MEB equation 
in Sec. III. The CMB and baryon ion-acoustic wave equations are then derived, 
respectively, in Sec. IV and V. Finally, we discuss the significance of plasma effects 
on cosmological perturbation. 

2. Plasma State of the Universe 

By saying " The universe is in plasma state" , we mean that relevant entities in 
the universe satisfy plasma conditions and hence the universe can be modelled 
as a plasma. As is said that a plasma is a quasineutral gas of charged and neu- 
tral particles which exhibits collective behavior^, two criteria-quasineutrality and 
collectiveness- must be fulfilled for the universe. As a perfectly uniform plasma is 
charge neutral and has no net electrostatic potential due to the perfect Debye shield- 
ing, the inevitable thermal fluctuations of its electron density at finite temperature 
would, however, induce a non-vanishing residual electrostatic potential of the ion 
that defies the Debye shielding. Such a residual electric force exerts an additional 
pressure on ions (baryons) which will be incorperated into the CPT to display the 
ion-acoustic behavior in the MEB equation in the following sections. The quasineu- 
trality now means that the scale of interests (i.e. the size of the universe) is much 
larger than than the Debye length: 



A- l/fco. io28(l + 2)-icm> Ad = VT^fceiyi^re^^T^. (1) 
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This makes the system is neutral enough so that one can take rii = n,, and not so 
neutral that all interesting EM interactions vanish. In addition, collective behaviors 
require enough particles inside the Debye sphere: 

Nd = n--KX% > 1- (2) 

Before decoupling {z > 1000), T^ = T^ ^ r^,o(l + z) and rig - IQ-^^n^ = 
10~"'^''n^^o(l + zY with the present photon temperature and number density T^.o = 
2.TiK and n^o — 400cm~'^ and redshift z > 1000. These information renders 

Ad = 1.5 X 10'*(l + z)^^ cm, (3) 

Nd = 5.7 X 10^ cm"''' (4) 

After reionization, T^ ^ Ti, — Tj^jl + z)^/(l + z* )^. As indica ted observationally, 
the redshift of reionization z* ranges as 5 < z* < 20^ I I I land the reionization 
temperature T^* is at least higher than lOOOiiT. Similarly, 

Xd = 1.6 X W^T^il + z)"^/^(l + z*)-^ cm, (5) 

Nd ^ 6.8 X 10^r(,*(l + zf^^{l + z*f cm^^ (6) 

The plasma conditions are apparently satified in both pre-decoupling and post- 
reionization era and, therefore, our plasma treatment is applicable. 

3. The Maxwell-Einstein-Boltzmann Equation 

Before decoupling, the EB equations can be derived systematically ■ - to account 
for the couplings between cosmological components. The equation displaying the 
acoustic feature on CMB spectrum is then rendered from these coupled equation. 
Following this systematical way, our MEB equation is deduced by introducing the 
EM force appropriately and combining it with gravity for the Boltzmann equation 
of baryons. 

The most general form of Boltzmann equation is 

'i-C[fl (7) 

where / is the distribution function and C'[f] characterizes the collision terms. In 
CPT, the distribution function / can be expressed as f{t,E{p),p^,x'^) with energy 
E and momentum p, moving in direction p at {x,t). For radiations (relativistic 
species), E ~ p and E^ ~ p^ + m^ for matter with mass m. Since the collective 
behavior emerges in the EB equation for baryons and photon EB equations remain 
unchanged, we present only the derivation of modification for baryons. Following 
the convention as in Ref. [22l the LHS of the Boltzmann equation is written as 

di^di didx^ didE didi_ 

dt dt dx' dt dE dt dpi dt ' ^ ' 

In equilibrium, the distribution function depends only on the length of the momen- 
tum, i.e. the energy E{p) and not on its direction p', so its dependence on p^ should 
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emerge from pertm'bations and be of the order. Therefore df /dp^ and dp^/dt are 
both first oorder terms and the last term in Eq. ([8]) is negligible in the linear theory. 
Choosing the conformal Newtonian Gauge, ds^ = —(1 + 2'^)dt^ + a^(l + 
2^)5ijdx^dx^ , where a is the scale factor, ^ the Newtonian gravity and $ the space 
curvature, the constraint gfi^P^P'^ = — w^ renders the four-momentum for baryons 

1 - $ 

E{l-^),pp' 



(9) 



where P^ = dx>^ /dS,, with S, parametrizing the particle's path. 
This definition renders directly 

^ = ^ = ^ = ^(1-* + $) (10) 

dE/dt (or dp/dt) is deduced from the geodesic equation 

^ = -r::, — — or,eqmvalently, ^ ^ -r::,P"P'^. (11) 

This geodesic can be regarded as a variation of Newton's second law describing the 
motion under the gravity in the form of metric perturbation. Since baryons feel not 
only gravity but also EM interaction, the geodesic equation should be modified to 
address the EM field as well. The equation of motion under Lorentz' force reads^ 

«'^-^--^"^^--^"^^"§' (12) 

d(;^ rrib dQ rrii, dt, 

where q is the charge of the baryon(proton) and P^^ = 9^A^ — di^A^^ is the EM 

tensor with four-potential A^. 

Here C, the proper time, is related to ^ via the constraint of P^ as 



d^ d^ d^2 ^^2- 



mi = P^P^ = -—-^ = — = --^. (13) 



Therefore d(^ = mbd^. 

By superposition principle, the geodesic equation is modified to address both gravity 

and EM force 

r/P'' 

-^^-KpP"P^ + lFSP"- (14) 

As the Christoffel symbol F^ „ is expressed in terms of g^^ 

^»f3 = i^g^'^idagup + dpga^ - d^gap), (15) 

the EM tensor P^i^, in terms of A^^, is derived from the electromagnetic lagrangian 

^EM - -^F'^-'F^, = 1^(^' - ^')' (16) 

OTT 47r 

where £ and B arc electric and magnetic fields induced from the charge density 
fluctuation, and related with electric and magnetic potentials and Ahy £ — —Vcf) 
and B = y X A, with the gradient operator Vi = ^^^ . 
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Evaluating Eq. p6|) renders the four-potential 



-(!-* + $), 

a a 



(17) 



As /x = designated and d^ replaced by dt, expanding Eg. p^ renders, with Eg. pU)) . 



di^di dj_p^_dj_l^p^ d^p^ pp^d^ 

dt dt dx^ aE dE\ E dt E a dx 



^ + 'S^a? • '"" 



the LHS of Eq.©. 

The last term inside the bracket describes the new effect from perturbed electric 
fields ensured by the total neutrality of the universe while all others are conventional 
terms. 0, like ^, is of the first order itself since only its gradient counts. The absence 
of vector potentials results straightly from the fact that they change the momentum 
but the energy of charged particles. Because of the p^ dependence, the extra em-term 
makes no contribution in the zeroth moment equation, like the curvature term. It 
contributes in the first moment equation and renders the MEB equation for baryons 



Vb + aHvb + iA;(* + -^) = — (SiOi + Vb), 

TTlb R 



(19) 



where "•" denotes derivatives with respect to the conformal time and R = |^ is 
the baryon-photon energy density ratio, 0i the first moment of the temperature 
perturbation O, and r the optical depth. 

Since the perturbed electric field 4> arises from the density perturbation, 



e(f) 


1 Awe^no ^ 

TjO — 


Ofc = 


ksT 


_ 


- 


nib 


nib k^ + A^ 


TTlb 


nib 



360, 



(20) 



where the last equality is enforced by the adiabatic condition. The MEB equation 



Eg. pg)) is then written as 

Vb — -3iQi + — 

T 



Vb + aHvb + zfc * 



knT 



SGr 



(21) 



4. CMB Anisotropy Before Decoupling 

The photon distribution function, with the temperature anisotropy = 6T/T, is 

-1 



f{t,p,p\x) 



exp 



P 



1 



(22) 



Tit)[l + e{t,p\x)] 
Since photons are not disturbed by EM interactions, the photon MEB equation 



is the same as its EB version 

e + ik^e + $ + i/c/i* = -r[eo - e + ^u^], 

where the moments are defined as, with Legendre polynomials P;, 

1 f^ dfi 



e, 



HY J-i 2 



-pKm)0(m)- 



(23) 

(24) 
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In the tightly coupled limit, the only nonnegligible moments are the monopole and 
dipole. Eq. p3|) renders the coupled equations of the zeroth and first moments, 






eo + fcei = -$; (25) 

ivb 



^^ 3 



(26) 



Eas. ([2Tl [25l and [26]) allow to extract our MEB equation for the temperature 
anisotropy 0o : 



Go + aH-^Qo + k" 
i + ri 



1 R fceT 



3(1 + i?) l + R m 



eo 



= 4v,-.H_A-<i,-^ (27) 

From this MEB equation, the sound speed of the CMB fluid is changed as 

^^^ : 1 I ^ ^bT 

' 3(1 + i?) 3(l + i?) 1 + i? mfc ' ^ ' 

where the former term represents the conventional theory and the last reveals the 
plasma collective contribution. The contribution to the sound speed from the EM 
induced pressure relative to the photon pressure is easily estimated 



r^RL ^ 10-5 (29) 

Accordingly, the sound horizon is changed to the same order. The plasma contri- 
bution should therefore occur at scales much larger than the Silk damping and is 
observationally negligible. 



5. Baryon Perturbation After Reionization 

After reionization, the universe is transferred to another plasma state and EM 
interactions influence the evolution of the baryon density perturbation again. In 
this epoch, we study the evolution of perturbations, using Newtonian gravity, as in 
Ref. [211 Since photons were decoupled long ago, Boltzmann equation in the integral 
form, the fluid equation is appropriate to describe the density pertubation. With 
EM forces included, our evolution equation for baryon perturbation is derived from 
plasma equations in the expanding universe 

^ = ^ + (ub.V)p, = -3i/p,; (30) 

^ = -VM'g, - -VP - -V0em, (31) 

at Pi, m 

where pb and Ub denote the mass density and velocity of the baryon and $gr, P and 
^om are the gravity, pressure and em-potentials the baryon suffers. The gravitational 
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potential should account for total mass pm including baryons and the dark matter 
via Poisson equation 

V^^gr = AnGp. (32) 

The universe is charge neutral in total so the electric field arises from the density 
fluctuation and is perturbative, (jjcm — 4'- The expansion is defined by the local 
Hubble parameter as 

H=^(V-Ub). (33) 

Perturbing and linearizing these equations render, to first order, 

{5pb)- = -ipbSH - SHSpb, (34) 

vt + Hvb = -\/^-—\7SP-—V(j}, (35) 

Pb mb 

6H=^\/- vb, (36) 

V^* = AnGSp, (37) 

where " • " denotes partial differentiation with respect to time t. Instead of Spb, the 
density contrast Sb = Spb/ Pb is more useful to work with. The MEB equation is 
then extracted from these coupled equations 

Sb + 2H6b - -\/^SP - — V20 = iirCpSb. (38) 

p nib 

Under Fourier transformation, our MEB equation for baryons in rcionization era is 

k,b + 2i?Jk,b + (-)'(l + h'^Si,,b - 47rGp(5k,;„ (39) 

Va/ 3 iTib 

where thermodynamics leads the SP-teim to 

5Pk Spi^fi SPk 5kBT 

= T = -li Kb- (40) 

Pb Pb opi^,b orub 

We therefore deliver the origin of our MEB equation, according to which our 
investigation of plasma suppression was performed in Ref. [IJ 

6. Conclusion 

The acoustic behavior is an impressive feature in the CMB spectrum and matter 
power spectrum. Our study shows manifestly the EM fields also contribute to this 
acoustic oscillation collectively. Though the whole system is neutral in total, the 
nonzero temperature induced charge density fiuctuations form a collective electric 
field which turns to be an extra portion of pressure exerting on the densithy per- 
turbations in both cases. However, the contribution on the CMB is too small to be 
observed. Both the order of magnitude estimation and the ion fluid equation analy- 
sis tell that the EM contribution is of the same order as of the pressure from charged 
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particles. As the adiabatic condition indicates baryon and photon density contrasts 
being parallel, the baryon to photon number ratio rih/n^ ~ 10^^*^ states that the 
EM effect on the sound speed is ~ 10~^ smaller than ordinary CMB fluctuation. 
Our generalization of the Boltzmann equation proves this number by incorporating 
the Lorentz force into the ordinary geodesic equation for baryons. 

Though the plasma collective effect is negligibly small in CMB spectrum, it 
should has sizable imprints in the matter power spectrum as discussed in our worleJ. 
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